Effects of nematic fluctuations on the elastic properties of iron arsenide 

superconductors 
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We demonstrate that the changes in the elastic properties of the FeAs systems, as seen in our 
resonant ultrasound spectroscopy data, can be naturally understood in terms of fluctuations of 
emerging nematic degrees of freedom. Both the softening of the lattice in the normal, tetragonal 
phase as well as its hardening in the superconducting phase are consistently described by our model. 
Our results confirm the view that structural order is induced by magnetic fluctuations. 

PACS numbers: 74.70.Xa ; 74.40.Kb ; 74.25.Ld 



In the recently discovered FeAs superconductors 0, Q , 
the rather large superconducting transition temperatures 
are in close proximity to an antiferromagnetic (AFM) 
phase transition at Tn and a structural transition from 
tetragonal (Tet) to orthorhombic (Ort) at T s [3-6]. This 
begs the question to what extent spin and lattice de- 
grees of freedom play important roles for high-T c super- 
conductivity (SC). The fact that Ts and Tn track each 
other closely (see inset of Fig. 3), as well as the exis- 
tence of an isotope effect for Tn [7|, are clear evidence 
that they signal the onset of strongly coupled ordered 
states. Since Tg > Tn, it seems, at first glance, plau- 
sible to assume that lattice effects are primary and the 
AFM order is simply a secondary effect. However, it has 
been suggested [8, 9] that spin fluctuations, related to the 
observed AFM ordering at Tn, will, at higher tempera- 
tures, lead to emergent nematic degrees of freedom (lp| 
that couple to the lattice 

[HG3. Even though of mag- 
netic origin, probing nematic fluctuations is nontrivial, 
since the magnetic field fluctuations associated to them 
average to zero. Thus, determining to what extent these 
fluctuations are present is key for identifying the relevant 
low energy modes in the pnictides. 

Recent experiments |l3l . found strong anisotropy in 
the electronic properties of the Ort phase, suggesting an 
underlying electronic nematic state. In this Letter, we fo- 
cus on the Tet phase and demonstrate that the measure- 
ment of the shear modulus C s allows for a direct deter- 
mination of the nematic order-parameter susceptibility 
Xneirf In particular, we show that the high-temperature 
shear modulus C Sj o is renormalized by nematic fluctua- 
tions to: 

C s 1 = C s q + A 2 C S q Xncm, (1) 

where A is the magneto-elastic coupling. We calculate 



Figure 1: (left panel) In-plane magnetic ordering of the pnic- 
tides. The two iron atoms present in one tetragonal unit 
cell (open blue circles and filled red circles) form two coupled 
AFM sublattices with co-planar staggered magnetization mi 
and rri2. (right panel) In the disordered phase, the spins in 
the two sublattices fluctuate around two arbitrary staggered 
magnetization directions. In the nematic phase, these two 
Neel vectors are locked either parallel (tp > 0) or anti-parallel 
(tp < 0) to each other, but still with (mi) = ( m 2) = 0. The 
coupling between the nematic and the elastic order parame- 
ters makes the bonds between neighboring parallel (antipar- 
allel) spins contract (expand). 



X ncm in the Tet phase using a 1/N approach and com- 
pare the results with Resonant Ultrasound Spectroscopy 
(RUS), which measures C s . Our RUS data show that 
nematic fluctuations are prominent over a large portion 
of the phase diagram of the pnictides. Although driven 
by AFM fluctuations, they represent new emergent col- 
lective degrees of freedom. Their coupling to the lat- 
tice makes the structural transition coincident with the 
nematic transition, implying that the elastic degrees of 
freedom are secondary with respect to the electronic ones. 
Here, we analyze the behavior of nematic fluctuations not 
only in the normal state, but also in the vicinity of T c , 
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explaining the observed hardening of C s in the SC state. 

In the AFM phase, spins are coupled antiferromagnet- 
ically for Fe-Fe bonds along one diagonal and ferromag- 
netically along the other diagonal of the two-Fe unit cell 
(Fig. 1). The in-plane magnetic ordering can be divided 
into two sublattices [lH with Neel magnetizations mi 
and rri2, which are weakly coupled since in the classical 
ground state, local fields associated with one sublattice 
cancel in the other. 

This sublattices coupling is related to a Z2 (i.e. dis- 
crete, Ising type) symmetry: mi — > —mi, 1112 — > m.2 
along with a rotation of the coordinates by ir/2. The 
Z2 symmetry, manifest by the two possible orientations 
of the magnetic stripes, will give rise to an Ising-like or- 
der parameter, which in the continuous limit is given by 
if (x) = mi (x) ■ ni2 (x) [10], with x = (r,r) as shown in 
Fig. 1. It describes the relative orientation of mi and 
1112 and does not change sign upon magnetic field inver- 
sion. The particular magnetic structure of the pnictides 
naturally accounts for this emergent degree of freedom, 
regardless of whether one considers an itinerant TfJ 17 1 or 
a localized [18| picture. Once (92) ^ 0, the Z2 symmetry 
is spontaneously broken [10], and so is the underlying C4 
symmetry of the lattice. This is why one associates the 
Ising-ordered phase to a "nematic" state, even though, 
strictly speaking, the full rotational symmetry, as well as 
the translational symmetry, are already broken above T$. 

To probe ip, we use the fact that the nematic order cou- 
ples to the elastic degrees of freedom of the lattice [lll[l9| : 
in orthorhombic fluctuations, present for T > Tg, bonds 
connecting antiparallel spins are ex pan ded whereas those 
connecting parallel spins contract |H El (see Fig. 1). 
There is therefore an energy cost associated with such 
an orthorhombic fluctuation, characterized by the shear 
modulus C s = Cqq. Shear strain (e s = 2e xy , where e a p is 
the strain tensor) and the nematic order parameter are 
linearly coupled in the energy: 



^el— mag 



= A 



(x) if (x) . 



(2) 



i.e. they order simultaneously. Here A > is 
the magneto-elastic coupling constant and J ... = 

rp-l 

J Q dr J d 3 r.... Assuming a harmonic lattice for A = 0, 
with bare (high-temperature) shear modulus C s ,o, we 
obtain the renormalized C s in Eq. (jTJ) . Note that 
the static nematic susceptibility Xncm = Xncm {q 0) 
is a four-spin correlation function, given by Xncm (q) = 
(ipqtf—q) — (ip q ) , where ip q denotes the Fourier transform 
of ip(x). From Eq[IJ we see that nematic fluctuations 
soften the lattice in the Tet phase (T > Tg). A diver- 
gence of Xncm brings C s to zero, signaling the onset of 
the Tet-Ort transition. 

To probe the shear modulus in the pnictides, RUS 
[22| was performed on undoped BaFe2As2 and optimally 
doped BaFei.84Coo.i6As2 single crystals grown out of 



FeAs flux [23]. More than 20 mechanical resonance fre- 
quencies / were detected. Since the crystals have 6 inde- 
pendent elastic constants C%j in the Tet phase, it is clear 
that each resonance line / will be proportional to a com- 
bination of Cij . Extraction of the full elastic tensor is not 
possible due to samples size limitations, but information 
about C s — Cqq can be obtained. Because the samples 
undergo a Tet-Ort transition, C s must get soft by sym- 
metry. Numerical simulations of the resonance spectrum 
show that this soft C, dominates several of the measured 



lines, with f 2 oc C s [22 1 



Following this procedure we present (Fig. 2) the T de- 
pendence of two resonance frequencies that are represen- 
tative of a larger set of C s -dominated lines. The observed 
reduction in f 2 indicates a dramatic softening of C s . In 
the undoped case (Fig. 2a), this softening is cut off by 
a simultaneous structural/ AFM weak first order transi- 
tion at T$ = 130K. In the optimally doped case (Fig. 
2b), where no AFM or structural transitions are present, 
the softening continues to lower temperatures until it is 
truncated by the SC transition at T c = 22K. In order to 
confirm that these features are not due to any particular 
mechanical instability of the samples, we also performed 
RUS in strongly doped samples that present neither SC 
nor structural transitions. In this case, only an uniform, 
mild increase in / with decreasing T was observed. 

Next, we calculate the nematic susceptibility Xnem in 
Eqffl Following Refs. @, [§] we introduce an action 
S = S'mag + S c i + S'ei-mag for the magnetic and elas- 
tic degrees of freedom, with S'ei-mag given by EqJ5J For 
each sublattice in Fig. 1, we consider: 



^mag [ m i\ 



1 



J X 



(3) 



where w > is a coupling constant and q = (q,iu) n ) refers 
to momentum, q, measured relative to the ordering vec- 
tor Q = (-,±— ), and bosonic Matsubara frequency u> n , 

with J q ... = TY^ n J The bare spin susceptibil- 

ity Xo,q i s assumed to be %o \ = r o + qf — 2 J z cos (q z c) + 
II(w n ). Here, the particle-hole bubble H(uj n ) is deter- 
mined by the electronic structure and leads to a sensitiv- 
ity of xo,q with respect to SC. Not only the static part of 
Xo,q is changed in the SC state, but also its dynamics j24[: 
the opening of the SC gap A makes the spin dynamics 
ballistic for u„<A(z=l), while for w„ ^> A one recov- 
ers the overdamped dynamics of a metallic AFM (2 = 2). 

Thus, we can write n (w n ) — n (0) = T z \u> n \ 2 ^ z with Lan- 
dau particle-hole damping coefficient Y z= 2 — 7 -1 as well 
as r z= i ~ (A7) -1 . ro characterizes the distance to a 
magnetic critical point and J z is a magnetic coupling be- 
tween Fe-As layers separated by c. To describe systems 
of different Co doping, we vary 7 and ro- The full mag- 
netic action is given by the sum of the contributions of 
the two sublattices and the dominant [j| [IlJ term that 
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Figure 2: Temperature dependence of the RUS measured 
squared resonant frequency f 2 (red points), and of the calcu- 
lated shear modulus C a = Cm of the tetragonal phase (solid 
blue line), for both (a) undoped BaFe2As2 and (b) optimally 
doped BaFei.84Coo.i6As2 . 



couples them: 



= ^5i°) g [m i ]-|^(m 1 (z).m 2 



(z)) 2 (4) 



where g > is the sublattice coupling. In a localized 
picture, q is related to quantum and thermal spin fluctu- 
ations [10( , whereas in an itinerant picture it is related to 
the ellipticity of the electron band 16, 25|. To solve for 
the coupled magnetic and elastic system we assume that 
rrii is an A^-component vector and analyze the leading 
term for large- N We find: 



Xnom (q) 



(<?) 



1 



9rX0, 



(?) 



(5) 



where g r = g + X 2 /C s ,o is the sublattice coupling renor- 
malized by magneto-elastic interactions. Thus, even if 
g = 0, the magneto-elastic coupling gives an effective cou- 
pling between the sublattices. Xo,nem (?) = N J X P Xv+q 
is determined by the dynamic magnetic susceptibility \ p 
of EqEJbut with tq +11 (0) replaced by r — . £ is the 
magnetic correlation length, determined in the large N 
approach. If d + z < 4 it follows that xo.ncm (0) diverges 



as £ — > oo. Here, the system has d = 2 (d = 3) behav- 
ior for r ^> J z (r -C Jz)- Thus, a sufficiently large, but 
finite £ suffices to cause a divergence of Xncm, implying 
that the Tet-Ort transition is induced by AFM fluctua- 
tions, not AFM order; hence Tg > 7n (except when both 
transitions are first order). 

First, we present our results for the Tet and non-SC 
state in Fig. 2, together with the RUS data. The spin 
dynamics is overdamped (Le. z — 2), as seen in in- 
elastic neutron scattering [26j. We used the parameters 
Tundoped = 2-2 x 10- 3 , tJoU = ^ x 10" 3 , g - 0.015e , 

— ~ = 10 _3 eo, and J z = 10 _3 eo, where e ~ 3.4 meV is a 
magnetic energy scale. The value for the Landau damp- 
ing is consistent with that used in Ref. [26j]. The lattice 
softening, seen in the experiment over a large tempera- 
ture range, can be well described by the above theory, 
suggesting nematic fluctuations up to room temperature 
and over a significant portion of the phase diagram. The 
amplitude of the effect indicates a strong magneto-elastic 
coupling, consistent with the observed softening of some 
phonon frequencies [27|. 

For both samples, the lattice softening in the RUS data 
stops at low T. For zero doping this happens at the joint 
first order structural/ AFM transition, where a hardening 
of the lattice is expected. Due to the Fermi surface re- 
construction in the AFM phase, however, a quantitative 
description of C s below T/v is beyond the scope of our 
model. 

For optimal doping, the lattice also gets harder below 
T c . In this case, though, there is a priori no reason for 
C s to be minimum at T c . Note that this behavior is 
distinct from the Testardi-behavior seen in other super- 
conductors, where strain couples to the squared SC or- 
der parameter, l 1 ^ 2 , and changes the transition temper- 
ature T c — > T c (s a b), making the elastic modulus behave 
very similar to minus the heat capacity at the transition 
28l 29| . Here, symmetry prohibits such a coupling and 
non-Testardi behavior is expected at T c . 

To understand our results close to T c we determine 
the dynamic spin response in the normal and SC phase. 
Analyzing the spin dynamics, ballistic for u <C A and 
overdamped for uj 3> A, we find that the frequency in- 
tegral in Xo,nem is dominated by overdamped dynamics 
(z = 2). Therefore, the key impact of the SC phase is 
not on the dynamics, but on the static part of the sus- 
ceptibility. Magnetism in the pnictides has been shown 
to be strongly affected by SC, to the extent that the 
AFM order parameter is suppressed below T c (30l . l3lj |. 
A Landau expansion for the competition of the AFM 
and SC order parameter is governed by the coupling 
Y I* (x) | 2 (m 2 (x) + ml (x)) [H [33j|. Close to T c this 
coupling leads to an increase in the inverse magnetic 

correlation length r — > r + A m ^l^ 7 (x)| 2 ^, which sup- 
presses the static part of the spin susceptibility x p . Thus, 
magnetic order and correlations are weakened in the SC 
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dom are important low energy excitations in a large tem- 
perature and doping regime of the system. 
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Figure 3: Shear modulus C s in the tetragonal, SC phase. 
The inset shows a schematic (x,T) phase diagram for 
Ba(Fei_ a; Co :E )2As2, with the green arrow denoting optimal 
doping. In the main panel, the solid line refers to the opti- 
mally doped system, whereas the dashed and dotted curves 
refer to systems that are respectively deeper in the overdoped 
region. The red points refer to the RUS data of Fig. 2b. 

state and xo.ncm (0) = N J p xf, acquires a minimum at 
T c . Given the overdamped spin dynamics we obtain 

Xo.nom (0) oc — In + A m ^|\I/| 2 ^ for 2D magnetic fluc- 
tuations. To show that this effect explains the observed 
anomaly of C s at T c , we compare in Fig. 3 the RUS data 
for the optimally doped sample with our theory. In gen- 
eral, the correction to £ -2 is proportional to the compos- 
ite operator ^| V I'| 2 ^ oc (T c — T) 1 a , yet the SC transition 

is well described by a mean- field approach (a = 0). On 
general grounds, one expects the AFM fluctuations to get 
weaker along the overdoped region, which would lead to 
a less pronounced softening of the lattice for T > T c . 

We note that the reduction of C s close to the Tet-Ort 
transition is a general feature, expected in any critical 
theory. However, the hardening of the lattice below T c , 
where there is no structural instability, is an unexpected 
observation. Yet, it finds a natural explanation in terms 
of magnetically driven nematic fluctuations and compet- 
ing AFM and SC order. Furthermore, the same model 
is able to explain the reduction of orthorhombic (i.e. ne- 
matic) order in the SC state, as observed in Ref. [HI]. 
The existence of nematic fluctuations is crucial to our 
results: if one considered only generic magnetic fluctua- 
tions, then the relevant magneto-elastic coupling would 
be m 2 £ii, instead of Eq. [2j In this case, the magnetiza- 
tion would couple not to the shear distortion, but to the 
longitudinal modes, which cannot explain the strong cor- 
relation between the AFM and the Tet-Ort transitions. 

In summary, the spin ordering in the pnictides re- 
sults in emergent nematic degrees of freedom that can 
be probed through the shear modulus. Our results show 
that both magnetic and nematic/shear degrees of free- 
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